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Abstract: A two-site Ising model is implemented as an injection-locked 
laser network consisting of a single master laser and two mutually coupled 
slave lasers. We observed ferromagnetic and antiferromagnetic orders in 
the in-phase and out-of-phase couplings between the two slave lasers. Their 
phase difference is locked to either 0 or tt even if the coupling path is 
continuously modulated. The system automatically selects the oscillation 
frequency to satisfy the in-phase or out-of-phase coupling condition, when 
the mutual coupling dominates over the injection-locking by the master 
laser. 
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1. Introduction 

Combinational optimization problems are ubiquitous in our modern life. Classic examples 
include protein folding in life-science and medicine, frequency assignment in wireless com¬ 
munications, navigation and routing, microprocessor circuit design, computer vision and voice 
recognition in machine learning and control in social network. Those combinatorial optimiza¬ 
tion problems often belong to NP (Non-deterministic Polynomial-time) hard problems, for 
which modern digital computers and future quantum computers cannot find exact solutions 
in polynomial time. 

Three-dimensional Ising models and two-dimensional Ising models with a Zeeman term are 
also known to be NP complete/hard problems ifTl. Quantum annealing has been proposed to 
solve such Ising models by utilizing quantum mechanical tunneling induced by a transverse 
field ||2]- m. This concept was implemented with superconducting flux qubits in the D-Wave 
One / Two machines 0. The idea of quantum simulation is to use a well-controlled quan¬ 
tum system to simulate the behavior of another quantum system i). Recent experiments using 
trapped ions ITl- lfTOll have realized quantum simulation of Ising models. 

We have proposed an alternative scheme for implementing such Ising models with an 
injection-locked laser network CUO . Another physical implementation of Ising model based 
on degenerate optical parametric oscillators is also proposed ifTSll where the binary Ising spins 
are represented by the bistable phase (0 or n) of a degenerate optical parametric oscillator 
(DOPO) and experimentally demonstrated m with M = 4 DOPO pulses. In principle, 
such injection-locked laser network and parametric oscillation network can be operated at 
room temperature and constructed as relatively compact systems. The Ising spins in the laser 
network studied in the present work are represented by the discrete phases or — ^) of 
single-mode slave lasers with respect to the master laser phase rather than two orthogonal 
polarization modes proposed in the original paper im. Here, we describe the first experimental 
implementation of a two-site Ising machine using an injection-locked laser network. 
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Fig. 1. Ising spin representation by using polarization state (a) or phase state (b) of slave 
lasers 


2. Ising spins in injection-locked slave lasers 

Our goal is to implement the Ising models with a Zeeman term in an injection-locked laser 
network. The Ising Hamiltonian with a Zeeman term is described as 

M M 

jr = + 52 A,C7 ,j (1) 

'<;■ i 

In this Hamiltonian, describes an Ising spin, i.e., spin projection onto the z-axis. is the 
interaction coefficient between spin i and spin j, and Z, is a supplemental Zeeman (external 
field) term. In contrast to the standard Ising model with nearest neighbour coupling, here we 
deal with an arbitrary graph in which a vertex can be connected to any other vertices. The 
injection-locked laser network proposed in our previous paper mu [El can find the ground 
state of the Hamiltonian, Eq. ([T]), through a laser phase transition. A photon in the lasing mode 
is not localized in any specific slave laser but its wavefunction is coherently spread over all 
slave lasers as partial waves of each individual photon. At the end of the computation, the 
phase configuration of such partial waves is expected to represent the spin configuration of a 
particular ground state = 1 ^ M}. 

Ising spins are represented by the polarization states of the slave lasers as shown in 
Fig. [TJa), where \R) and \L) represent right circular (up-spin) and left circular (down-spin) 
polarizations, respectively. Each slave laser is injection locked by a master laser and prepared 
in the vertical polarized state |y). Then, the polarization state evolves towards either the right 
or left circular polarization after the mutual couplings between the slave lasers are turned on. 
Such polarization evolutions are described as a change in the relative phase between the equal 
amplitude diagonal polarization states \D) and \D) (see Eig. [TJa)). The amplitudes of two 
diagonal polarizations \D) and \D) are held constant but their phase difference continuously 
evolves during the computation. This phase-based picture, immediately makes it apparent that 
we need only one polarization state |Z)) because the phase of the another polarization state 
\D) is inversely-symmetric to the phase of |Z)) state. In the new model, the Ising spins are 
represented by the phase of the single slave laser mode, as shown in Eig. [TJb). The zero phase, 
as the reference, is determined by the phase of the injected master laser. When the master laser 
injection is reasonably strong and the self-oscillation frequency of the slave laser is identical 
to that of the master laser, the phase of the slave laser is identical to that of the master laser, 
i.e. 01, = 0. The phase rotates to 0^, = ±;r/2 from the initial phase 0^, = 0 due to the mutual 
coupling between the slave lasers , which corresponds to the Ising spin rotating to the up or 
down state. 







3. Experimental setup 
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Fig. 2. Experimental setup. The blue line shows the mutual coupling between two slave 
lasers (DMLDs). The red line shows uni-directional injection from the master laser that 
implements the Zeeman term The relative phase of DMLDl and 2 is measured by the 
Mach-Zehnder interferometer on the yellow lines with the FPI and APD. PDl, PD2, APD 
and FPI each have a refractive neutral density filter (~ lOdB attenuators) inserted in front 
of them as isolators. [Inset(a)]: How to implement the Ising interaction between two slave 
lasers. In-phase coupling between two slave lasers corresponds to ferromagnetic coupling. 
[Inset(b)]:Schematic picture of two site Ising model implementation. 

In our experimental setup, a master laser signal is uni-directionally injected into two 
slave lasers, which also implements the Zeeman term in Eq. 0 if the phase is slightly rotated 
towards either -F| ot from the standard phase zero. The slave lasers are mutually coupled 
by bi-directional injection-locking, which implements of the Ising coupling term in Eq. 0. 
The detailed experimental setup is depicted in Eig.|^ The optical path length between the two 
slave lasers and those between each slave and master laser are independently controlled by 
a polarization-dependent optical path. We used a tunable external cavity single-mode diode 
laser as the master laser(Koshin Kogaku LS-601A-56S2); it oscillates around 1577.5 nm with 
a spectral linewidth < 30 kHz. We used discrete mode diode lasers (DMLD, Eblana photonics 
ELDs) with a spectral linewidth < 100 kHz as the slave lasers. 

Two slave lasers are mutually coupled along the blue line in Eig. These slave lasers are 
nearly identical, and their injection currents are set to be the same in the experiment. How¬ 
ever, their temperatures must be tuned so as to make their self-oscillation frequencies identical. 
The optical path length between the slave lasers is L = 1550 mm. Moreover, their optical path 


























































length is controlled by piezoelectric transducer 1 (PZTl), wherein the modulation voltage is 
used to shift the position of the reflection mirror with a quarter wave plate 1 (QWPl). The role 
of the quarter wave plate 1 inserted before PZTl is to transform the vertical polarization light 
from the slave laser 1 to the horizontal polarization light before transmitting through PBS4. 
Likewise, the quarter wave plate 2 transforms the incident horizontal polarization light to the 
vertical polarization light before reflecting from PBS4. The same transformation is performed 
onto the output light from the slave laser 2. The two slave lasers oscillate in the horizontal polar¬ 
ization. The polarization beam splitter 1 (PBSl) and half wave plate 1 (HWPl) work together 
as attenuators for the beam from one slave laser (DMLD2), while PBS2 and HWP2 attenuate 
the light from the other slave laser (DMLDl). HWP3 is set at 0 = 45° to the polarization axis 
so that the polarization of the incident beam from the slave laser 1 (DMLDl) into HWP3 is 
rotated to the horizontal polarization and it is reflected to QWPl by PBS4. The master laser 
signal is injected into both slave lasers. The master laser oscillates in the vertical polarization 
and is converted into the horizontal polarization when it goes through HWP3. The optical path 
length for the master laser signal is not influenced by PZTL The phase of the master laser sig¬ 
nal injected into slave laser 1 which comes from the upper branch of the red line is modulated 
by the mirror mounted on PZT2. The output beams from the slave laser 1 and slave laser 2 
are partially reflected at BSl and BS2 and combined at BS3 along the yellow line. The PZT3 
modulates only optical path length from the slave laser to the BS3 without toughing the optical 
path length from the slave laser 1 to BS3, so that the interference pattern between the two slave 
laser outputs can be obtained at APD. The spectral linewidth and the interferometer output be¬ 
tween the two slave lasers are detected with a scanning Fabry-Perot interferometer (FPI) and 
an InGaAs/InP avalanche photo diode (APD). As shown in inset(a) of Fig|^, the ferromagnetic 
and the anti-ferromagnetic couplings are implemented by in-phase and out-of phase coupling 
path length between the two slave lasers.The inset(b) of Fig|^ shows a simplified picture of 
the experimental setup in Fig|^ A master laser signal is injected to the two slave lasers and 
a path from the master laser to the slave laser 1 is modulated by PZT2. The mutual coupling 
path between the two slave lasers is modulated by PZTL A path difference in the interference 
between the two slave lasers is modulated by PZT3. 


4. Slave laser phase modulation by injection-locking 

The phase of the injection-locked slave laser is continuously shifted from that of the master 
laser if there is a frequency detuning between the master laser and slave laser within the locking 
bandwidth m. The phase of an injection-locked slave laser shifts continuously from — f to ^ 
inside the locking bandwidth. Denoting the frequency detuning between the slave laser and the 
master laser by A(0, the relative phase shift (po satisfies the following relation: 
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where cOrQ is the self-oscillation frequency of the slave laser without injection-locking and co is 
the master laser oscillation frequency, Fq and Aq are, respectively, the injection held amplitude 
from the master laser and slave laser internal held amplitude, and Qg is an external cavity 
quality factor of the slave laser, a is the linewidth enhancement factor, j/pp is normalized to 
express a photon flux, while |Aop is normalized to represent a photon number. 

We measured the phase shift of the slave laser due to the frequency detuning A/ = Acojln 
between the master laser and the slave laser within the locking bandwidth. The frequency of 
the master laser is externally modulated at a 10-Hz repetition rate through its injection current. 
In Fig.j^the red line shows the applied voltage used to modulate the injection current. The blue 


line shows the interference signal between slave laser 1 and the master laser observed at BS5 
in Fig.|^for identical output powers of 1.12 mW. The range of the frequency excursion of the 
master laser is 4.5 GHz to cover the full locking bandwidth. The observed locking bandwidth 
is ~ 1.3 GHz (the locking-bandwidth boundaries are indicated by the green dotted lines). Since 
the slave laser has a negligible linewidth enhancement factor, a ~ 0, the locking bandwidth 

is given by A/lb = \/^ lUll H?) where = 7ja)|Tbp and Pout = are the 

input and output powers of the slave laser. However, in the slave laser (DMLD) used in this 
experiment, most of its cavity is formed by a passive waveguide without carriers so that the 
effective linewidth enhancement factor is negligible. This is experimentally manifested by the 
symmetric phase excursion of the slave laser for upward and downward frequency detuning 
shown in Fig. The phase of the slave laser is supposed to shift from to with the 
frequency detuning between the master and slave lasers El. The relative phase shift shown in 
Fig .[^is monotonical with the frequency detuning, which indicates the linewidth enhancement 
factor a of the slave laser (DMLD) is indeed small. 
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Fig. 3. The blue trace represents the Interference output between the master laser and slave 
laser 1 as a function of the frequency shift of the master laser. The red trace is the ramp 
voltage to sweep the master laser frequency. 


5. Numerical simulation for two mutually coupled slave lasers 

In the phase model of an injection-locked laser network, the binary phase of each slave laser 
acts as an Ising spin El- The c-number Langevin equations of the time-dependent amplitude 
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Fig. 4. Bifurcation of two slave lasers’ phases in mutually coupled semiconductor lasers: 
(a)The Ising term is dominant when the mutual coupling between the slave lasers is rel¬ 
atively strong (77 = 0.04) (b)The mutual coupling ratio for the Ising term is weaker than 
the case of (a) (rj = 0.01). (c)The master injection term is dominant when the mutual cou¬ 
pling is relatively small (rj = 0.005). The other numerical parameters are Ai = X 2 = 0, 
q}/Q = 10 l 2 ( s-i), Tsp = 10-®(s), C = 0.005 and /3 = lO^^ 


A;(f), phase and carrier number Nci{t) for an i-th slave laser are described as. 
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Here F^i, F^. and Ffji are the Langevin noise terms for the amplitude, phase and carrier number. 
CO/Q is the cavity photon decay rate, ^ is the optical attenuation coefficient of the master laser 
injection signal, rj is the common attenuation coefficient of the Zeeman term and Ising term, Tsp 
is the spontaneous emission lifetime, P is the pump rate into the slave laser i, Ecvi = PNi/Xsp 
is the spontaneous emission rate for i-th slave laser. j3 is the fractional coupling efficiency of 
spontaneous emission into a lasing mode in. 

Figure 1^ shows numerical simulation results for the time-dependent phases of the two slave 
lasers, where the Ising coupling J 12 is -Fl (anti-ferromagnetic coupling). In this case, the steady 
state solution is an anti-ferromagnetic order, which means the two slave lasers have opposite 
phases. In the numerical simulation, the slave laser phases are initially set to the reference phase 
(zero), which is determined by the master laser injection signal. When the mutual coupling term 
is switched on, as shown in Fig.j^a), the slave laser phases bifurcate into +j and — On the 
other hand, when the injection from the master laser is relatively strong compared to the mutual 
coupling, as shown in Fig.j^c), the slave laser phases stay in the reference phase (zero). When 
the mutual coupling between two slave lasers is slightly decreased, the relative phase difference 
between them is smaller than n as shown in Fig. I^b). It is assumed that the frequency of the 
two slave lasers are locked to the master laser frequency (O in these numerical simulations. 




















6. Switching between ferromagnetic and anti-ferromagnetic phase orders 


We experimentally studied how the slave laser phases are modulated when the optical path 
length between the two slave lasers (DMLDs) is varied. Notice that the optical path length 
modulates the polarity of the Ising coupling term Jtj. We observed a transition from one regime, 
where the Ising coupling is dominant, to the other regime, where the master injection term is 
dominant, when we changed the ratio of the uni-directional injection from the master laser into 
the two slave lasers to the mutual coupling between the two slave lasers by rotating HWPl 
or HWP2 in Fig. Figure (a) shows the interference signal between the two slave lasers 
(DMLDl and 2) when there is relatively strong mutual coupling between them, with PinMi = 
l.SfiW and Pi„Si 2 — 33/rW. Here P,„Mi is the injection power from the master laser into a fiber 
coupler of the slave laser 1 and Pi„Si 2 is the injection power from the slave laser 2 into the 
fiber coupler of the slave laser 1. It is observed that the strongly coupled slave lasers prefer 
to oscillate in discrete phases. The relative phase difference between two slave lasers seems 
to be either 0 or n, as demonstrated by an interference signal in Fig.|^a). The frequencies of 
the two slave lasers now deviate from that of the master laser, as described in the next section. 
We observed that the slave lasers still oscillate at identical frequencies, which is true even after 
the injection master signal is turned off. The voltage applied to PZTl linearly increases as a 
function of time in Fig. The PZTl displacement is Ad / AV =12.23±2.5 x 10^^ m/V with 
respect to the applied voltage AV. The optical path modulation for mutual coupling causes the 
phase difference between the two slave lasers to have discrete jumps every AV ^ 3.6V. The 
high and low outputs correspond to the ferromagnetic and anti-ferromagnetic phase orders. 
When the master laser injection and the mutual coupling signals are set to be PinMi = 2.2/iW 
and PinSi 2 — 03liW, the interference signal doesn’t show any modulation against the mutual 
coupling path modulation (Fig.j^c)). When the two signal intensities, are set to be on the same 
order (PmS^ — 2.0/iW) as shown in Fig.|^b), the interference signal shows a transition between 
the mutual coupling dominant regime in Fig. [^a) and the master injection dominant regime in 
Fig.^c). 
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Fig. 5. Interference signal output to show the relative phase A(j)i.i 2 between two slave lasers 
(DMLDl and 2) when the mutually coupling optical path length is modulated: (a) Ising 
coupling term is dominant, (c) Master injection term is dominant, (b) Transition regime 
between (a) and (c). 























Fig. 6. Interference signals to show the relative phase A(j>ri 2 between the two slave lasers 
(DMLDl and DMLD2). The green lines show the case that the interference path length is 
modulated. The blue lines show the case that the injection phase from the master laser to 
the slave laser 1 is modulated, (a) The Ising coupling term is dominant under a weak master 
signal injection, (b) The master signal is dominant under a strong master signal injection. 


Figure shows the interference output when the two slave laser signals are combined at 
BS3 and detected by the APD as shown in Fig.[^ The blue line shows the interference signal 
when PZT2 is electrically driven linearly as a function of time, which causes the phase of the 
master laser signal injected into the slave laser 1 (DMLDl) to be modulated. The blue trace 
in Fig. I^a) shows the case of using the same experimental parameters of Pi„Mi = 1.5/iW and 
PinS^^ = 33mW as Fig. ga), where the mutual coupling is stronger than the uni-directional 
master signal injection. Here, the phase modulation of the injected master laser signal is not 
effective. The green trace in Figj^a) shows the interference signal when only PZT3 is linearly 
modulated as a function of time, which manifests the phase coherence between the two slave 
lasers. The blue trace in Fig.j^b) shows the interference signal with PZT2 modulation under a 
relatively strong master signal injection (PmMi = 147/rW) compared with the mutual injection 
power {PinSi 2 = 33/fW). Here, the interference signal has the same phase modulation of the 
master laser. It suggests that the oscillating phase of slave lasers are mainly determined by the 
injected master laser signal and there’s phase coherence between two slave lasers as shown by 
the green trace in Fig.|^b). 


7. Spontaneous frequency selection of mutually coupled slave lasers 


PZT voltage AV(V) 


0 -10 - 20 - 30 - 40 - 50 



Time (ms) 

Fig. 7. Frequency shift of a slave laser due to the mutual coupling path length modulation 
as a function of time under a strong mutual injection (blue line) and a weak master signal 
injection (green line). 

We observed a sawtooth variation in the slave laser frequency when the optical path length 
between the slave lasers is modulated, as shown by blue dots in Fig.|7] Here the mutual coupling 
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Fig. 8. Schematic explanation for the sawtooth frequency modulation and phase jump in 
two coupled slave lasers, which explain the experimental results shown in Fig. and Fig. 
simultaneously. 


is dominant over the master laser injection. The experimental results shown in Fig.j^a) indicate 
that the oscillation frequencies of the two slave lasers are synchronized in the strong mutual 
coupling regime. We measured the frequency shift of the slave laser 1 (DMLDl) with the FPI 
at the top right comer in Fig. , while the signal from the slave laser 2 (DMLD2) is blocked 
just before BS3. Figure |^plots the slave laser frequency shift against the mutual coupling path 
length modulation between the two slave lasers. When the mutual coupling is stronger than the 
master laser signal injection, a sawtooth frequency shift of ^ 30 MFIz is observed as shown 
by the blue line in Fig. In this case, the phases are pinned either in the same phase or in 
the opposite phase, as shown in Fig.j^a). On the other hand, the oscillation frequency of the 
slave laser is held constant, when the master laser injection is dominant, as shown by the green 
line in Fig. Here discrete jumps in the oscillation frequency are experimental artifacts due to 
the discretization error(resolution limit) of the FPI. There are = 21 sawtooth-frequency-shift 
cycles when the applied voltage is swept over AV = 78V, so that frequency shift occurs at every 
applied voltage of AV = 3.7V, which corresponds to a PZTl shift of AcIl = 456.3 ±92.8nm. 
This sawtooth-frequency-shift is considered to be the result of spontaneous selection of the 
synchronized oscillation frequencies of the two slave lasers to minimize the overall threshold 
pump current or to maximize the output power ifT^ . The favorable phase difference between 
the two slave lasers is either 0i — 02 = 0 or n so that the two slave lasers build up a standing 
wave between the edges of the two slave lasers. 

Figure is a schematic explanation of the spontaneous frequency optimization. When the 
optical path length between the two slave lasers is an integer multiple of the self-oscillation 
wavelength Am = 2;rc/a), the number of loops of the standing wave is even, so that the two 
slave lasers oscillate with the same phase. In the second figure of Fig. [^a), the oscillation 
wavelength become slightly elongated as the optical path length L is extended while the phase 
difference between the two slave lasers is still — 02 = 0. However, when L is shifted 
from the resonant condition in the hrst figure of Fig. |^a), the oscillation frequency returns to 
the self-oscillation frequency (0 and the number of the loops of the standing wave between the 
two facets is odd, so that the phase difference between the two slave lasers is 0i — 02 = 
We plot that description of oscillation frequency shift as a function of L in Fig. [^b). N is the 
number of the loops and corresponding phase differences 0i — 02, are written at the top. The 
light blue line shows the slave laser oscillation frequency shift according to the linear shift of L. 
This frequency shift for a path length variation AL= ^ corresponds to the free spectral range 
^Vfsr = s the optical cavity length L. 

The path length variation for a single frequency shift event in our experiment as shown in 
Fig .|^is AL = ^ because the path length between the two slave lasers changes twice when the 
position of the mirror just before QWP is changed. The expected position shift of PZTl for a 












single sawtooth frequency shift is AL = ^ = 394.4nm which is close to the measured position 
shift of PZTl Adi = 456.3 ±92.8nm. The optical path length between the two slave lasers in 
Fig. I^is L = 1550mm. The sawtooth modulation in Fig. [^is ^ 30 MHz, which is less than 30 
percent of the estimated frequency shift AvpsR = ^ = 96.8MHz. The difference between the 
experimentally observed frequency excursion and the free spectral range stems from the phase 
shifts associated with the reflection from the two slave laser facets. As indicated in Eq. (3), the 
slave laser phase shifts away from the injection signal phase when there is a detuning between 
the self-oscillation frequency (Oro and the actual oscillation frequency (O. 

8. Conclusion 



Fig. 9. Ising machine based on a higher harmonic fiber mode-locked laser with A — 1 delay 
lines and the deferential phase detection. 

We have implemented a two-site Ising model by using the two mutually coupled slave 
lasers that are simultaneously injection-locked by a single master laser. As shown in Fig.j^a), 
the two slave lasers feature the ferromagnetic and anti-ferromagnetic orders when the two lasers 
are coupled in-phase (Jij < 0) or out-of-phase (Jij > 0), which is the ground state of the Ising 
Hamiltonian. The transition from the Ising-term dominant regime to the master-injection-term 
dominant regime was observed. The phase of the slave laser is continuously synchronized with 
that of a master laser when the master laser injection is stronger than the mutual coupling 
between the slave lasers. On the other hand, when the mutual coupling is stronger than the 
master laser injection, the phases of the two slave lasers jump discretely and the oscillation 
frequencies continuously shift with the optical path length modulation between the two slave 
lasers. In that case, the mutually coupled slave lasers communicate with each other and oscillate 
at the optimum frequency, where a standing wave develops between them that results in either 
a ferromagnetic or anti-ferromagnetic phase order. We note that as related work, the dynamic 
and chaotic oscillation behaviors of mutually coupled semiconductor lasers have been studied 
numerically and experimentally flUEol. 

For practical implementation of a larger-scale laser network based Ising machine, we can 
use a higher harmonic fiber mode-locked laser as shown in Fig. In this implementation, 
N optical pulses are simultaneously produced in a single fiber ring cavity so that they have 















identical oscillation frequencies. In order to implement the Ising coupling term, JijOi^ai^, we 
can introduce — 1 optical delay lines with independent intensity and phase modulators. The 
pick-off signal from i-th pulse is properly modulated in its intensity and phase, and coupled to 
the j-th pulse with an appropriate delay time. This time division multiplexing scheme allows 
us to implement N{N— l)/2 mutual coupling paths by only — 1 delay lines. Finally, the 
computational results are recorded by the one-bit delay differential phase detection ifiTll . When 
implementing with a higher harmonic mode-locked fiber laser with IGHz repetition frequency 
and 20m fiber ring cavity, the total number of optical pulses is ^100 in this system. If we 
increase the clock frequency up to lOGHz and the fiber length up to 200m, we can implement 
N ~ 10,000 optical pulses. 
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